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ABSTRACT 

A momentum dependent quark mass may be incorporated into a quark model in a manner 
consistent with dynamically broken chiral symmetry. We use this to study the high Q^ be- 
havior of the vector, axialvector, scalar and pseudoscalar two-point correlation functions. 
Expanding the results to order l/Q^, we show the correspondence between the dynamical 
quark mass and the vacuum condensates which appear in the operator product expansion 
of QCD. We recover the correct leading logarithmic Q^ dependence of the various terms 
in the OPE, but we also find substantial subleading corrections which are numerically 
huge in a specific case. We conclude by using the vector minus axialvector correlator to 
estimate the 7r+ — n^ electromagnetic mass difference. 



Introduction 

For large momenta, the QCD running coupling is small and calculations may be carried 
out systematically by a perturbative expansion in the coupling. As the momentum scale 
is lowered, nonperturbative effects become significant and the coupling expansion is no 
longer useful. Consider specifically the physics of light quark flavors in the context of a 
two-point correlation function, n(g^), which we define to be dimensionless. From a purely 
phenomenological standpoint, the onset of nonperturbative effects can be parameterized 
by a series of correction terms. 0] 
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Throughout this work, the tilde is used as a reminder that perturbative physics has been 
excluded. Q^ = —q^ is the (Euclidean) momentum squared and the C„ depend at most 
logarithmically on Q^. When Q^ (which we always choose to be positive) is not too small, 
this series is dominated by the first few terms and can be used to perform calculations at 
lower Q"^ than is possible using only the perturbative coupling expansion. (P is actually 
an operator product expansion (OPE) and the C„ are known functions of various nonzero 
vacuum condensates. |1[] 

The chiral symmetry breaking vacuum condensate with the lowest mass dimension is 
(gg), and it (when multiplied by a current quark mass) will contribute to C2 in (|l|). There 
is one other condensate in C2, (agG^G'^^'^), but it does not break chiral symmetry and 
will not be of direct interest in what follows. Since there are no condensates with mass 
dimension two, Ci=0 in QCD. 

The existence of a nonzero (qq) implies an effective quark mass0, and this in turn 
modifies the QCD Ward-Takahashi (WT) identities for the couplings of vector and axi- 
alvector fields to quarks. The complete two-point correlator in QCD is shown in Fig. ^ 
where one vertex is a bare 7^ [75] but the other vertex and both quark propagators rep- 
resent fully-dressed nonperturbative quantities that must satisfy the WT identities. In 
what follows we will consider only the contributions of the effective quark mass to these 
nonperturbative quantities, and it will be shown that for both the vector and axialvector 
correlators, a "minimal" WT vertex correctly reproduces Ci and the (qq) term of C2. 
We obtain expressions for C3 in the vector and axialvector cases which have the known 
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leading logarithmic running, and we point out the possibility of substantial corrections in 
C3 due to terms that are formally suppressed by logarithms. 

The same diagram (Fig. ^ also represents the scalar and pseudoscalar two-point cor- 
relators in QCD, although there are no WT identities to constrain the corresponding 
vertices. If the full scalar or pseudoscalar vertex is approximated by a bare one, C2 is not 
generated in its correct form. Improved scalar and pseudoscalar vertices can be found by 
appealing to the gauged nonlocal constituent (GNC) quark model. [|] 

The GNC Lagrangian contains constituent quarks with momentum dependent masses 
and pseudoscalar mesons, constructed to model dynamically broken chiral symmetry. 
The couplings of vector and axialvector fields to quarks are precisely the "minimal" WT 
vertices mentioned above. Although the GNC Lagrangian originally included the scalar 
and pseudoscalar fields in a trivial way|^], the Lagrangian allows for a natural extension 
of these couplings in a manner analogous to its vector and axialvector couplings. For 
the case of the scalar correlator we verify that this produces the correct expressions for 
the chiral symmetry breaking pieces of C2. C3 also has the correct leading logarithmic 
running, but here we find enormous subleading correction terms as well. 

It is important to stress that we will use the GNC model only to determine the form 
of the nonperturbative propagators and vertices appearing in the general correlator of 
QCD in Fig. |1|. Our main goal is to study the relation between the 1/Q^ expansion of 
these correlators and the momentum dependent mass function. We will need to consider 
the low energy behavior of these correlators only when we treat the pion mass difference. 
At low energies we may apply the GNC model directly, since it has been shown to model 
low energy phenomenology rather well0m. Thus at low energy the correlators will be 
described by the GNC model diagrams of Fig. ^ and at high energies they will be described 
by Fig. [|. The GNC model includes the pseudo-Goldstone bosons of QCD, and virtual 
effects of these mesons are accounted for in the low energy contribution according to the 
standard chiral Lagrangian approach. Meson loops are naturally cutoff in the model at 
the point where these particles lose their particle-like nature. 

The vector minus axialvector two-point correlator is of special interest in our study 
since it would vanish if chiral symmetry was not broken. We will choose an explicit 
form for the effective quark mass which becomes the known form[0] at large momentum 
scales and which resembles the successful00 GNC ansatz at small scales. We can then 
calculate the vector minus axialvector correlator numerically at any momentum scale by 



matching the low energy GNC model (Fig. |^ plus meson loops) to our high energy model 
(Fig. H) at an intermediate scale. An integral over all momenta produces the 7r+ — 7r° 
electromagnetic mass difference]^. 

The general 1/Q^ expansion 

It is convenient to express the vector, axialvector, scalar and pseudoscalar two-point 
correlators in terms of dimensionless functions of Q^, 

z|dSe^^-^(0|T{\/;(x)\/j'(0)}|0) = {q,q.-q^g,.YW\Q^) (2) 

ijd'xe^'^-^{Q\i:{Ai{x)Al{m^) = (gM9.-gV)^''n'^'(Q')-?/.9/'°n'^'(Q') (3) 

z/'d^a;e*«-^'(0|T{5"(a;)5^(0)}|0) = Q^ ^W\Q^) (4) 

Jd^xe'^-"(0|T{P"(x)P'(0)}|0) = g2Pn'^6(g2^ (5) 

where a, h are SU(A^/) flavor indices. Throughout this work we will restrict the discussion 
to light quarks. According to Fig. |l| we must determine the full quark propagator and the 
full vertex for each correlator. 

The most general quark propagator is 

If we think of writing the nonperturbative contributions to the inverse propagator as a 
1/Q^ series, analogous to (|T]), there are only two condensates with mass dimension less 
than five: (gg) and iasG'' 0"^^^"^ In the limit of vanishing current quark masses, we 
see on purely dimensional grounds that Z^Q"^) cannot contain (gg) while Ti{Q'^)/Z{Q'^) is 
independent of (osG^^^G'^'^'^). For our analysis of chiral symmetry breaking, we will set 
Z(Q^)=1 at the expense of omitting the gluon condensate (which contains no information 
about chiral symmetry breaking) as well as other effects at higher order in the 1/Q^ series. 
Recall that the vr"*" — tc^ electromagnetic mass difference is independent of (asG"^ G^' 
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We will retain S((5^) with the correct asymptotic behavior] 



S(g^) - ^^^ m , as Q^ - oo (7) 
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d = 12/(33 — 2Nf) for Nf quark flavors, /x is a renormalization scale, A is the QCD scale 
and the current quark masses are set to zero. With this normalization, (qq) represents a 
single quark flavor summed over three colors and summed over Dirac indices. 

The vector and axialvector vertices ^'^r^(p, q) are constrained to satisfy the Ward- 
Takahashi identities. 



q^''T';{p,p + q) 



S-\p + q)Y-Y^'\p) 
'5-'(p + g)75Y+75y5"'(p) 
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The A° are flavor generators normalized by Tr(A"A'') = 26"'^. These identities, plus the 
requirement that the vertices contain no unphysical singularities, uniquely define the 
longitudinal part of the vector and axialvector vertices. We will choose the minimal 
vertices by ignoring any extra transverse pieces that may exist. [^ The resulting vertices, 
for incoming(outgoing) quark momentum p{p' = p + q), are 



r;:(p,p') = y 
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r^(p,p') = V 
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Again, P^{P' ) = —p'^{—p' ). Notice that the vector vertex is completely free of sin- 
gularities (assuming none are contained within S(Q2)) but that the axialvector vertex 
is required to have a singular point. This massless state is the Goldstone boson of the 
dynamically broken symmetry — the pion for Nf=2. 

There are no analogous identities to constrain the forms of the scalar and pseudoscalar 
vertices. One might be tempted to adopt the bare vertices 



Sj^a 
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r"(p,p + g) = i—i^ 
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but this will lead to a disagreement with the OPE. One must also decide how to include 
current quark mass effects. Our approach will be to appeal to the GNC quark model 
which is the minimal Lagrangian that contains constituent quarks with mass S((5^) and 
which respects the dynamically broken chiral symmetry. The original GNC model was 
used at low energy scales, and the external fields were coupled to quarks according to 



( P^ - |T5|) , but as stated above, the scalar and pseudoscalar couplings contradict the OPE 
at larger scales. We therefore propose the following " GNC " model, which is identical 
to the original version except that the S and P fields now appear in the path-ordered 
exponential as well as in the local term. In Euclidean spacetime, 

— I i 

>CGNC'(a;, y) = ij{x)6{x - y)i^[d^ - iV^{y) - ii^A^ijj) + -i,,S{y) + -l^,l->P{y)]^{y) 

+V^(x)s(x - y)e(x)x(x, y)ay)i^{y) (i6) 
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X{x,y) = Pexp 



(17) 



'-e{zM - ^y,i^) + ^lrA,{z) - ll,S{z) + '-^^^,P{z)]^{z) (18) 



^(x) = exp 



/ a=l ^ 



(19) 



n"" contains the Nj-l pseudoscalar mesons with decay constant /, and ip is the A^j-plet 
of quark fields with mass T,{Q'^), the Fourier transform of S(x — y). V^, A^, S, P are 
the external fields. X{x,y) is a path-ordered exponential. Notice that we have assumed 
the same mass function S((5^) for all quark flavors, which means that (qq) is also flavor- 
independent. 

The Feynman rules for the GNC' model can be obtained from the Lagrangian (|16D in 
a systematic manner as described in ||^. Notice that the meson pole which contributes to 
the full axialvector and pseudoscalar vertices now appears explicitly, as shown in Fig. ^ 
Although the GNC Lagrangian does not contain an explicit meson propagator, a propa- 
gator is generated by the quark loops of Fig. H. The explicit calculation shows that the 



GNC Lagrangian does indeed reproduce the vector and axialvector vertices of (12-13). To 



include the effects of a nonzero current quark mass matrix, Ai, we make the substitution 

S{x) ^M + S{x) (20) 



For our purposes, we will only require terms linear in A4. 

With Feynman rules in hand, the expressions for the two-point correlators of Fig. |l] 
are easily written down in the form of 4-momentum integrals. The integrands are largest 
when the momentum-squared flowing through one quark propagator is of order S^(0) 
so that the other propagator has a momentum-squared of order Q^. When the results 
are expanded in powers of Ai and l/Q^, we obtain a series of the form (|l]). The simple 
procedure of expanding the integrand in 1/Q^ generates integral expressions for the various 
Cn, and these integrals become more divergent for increasing n. We stress that the full 
expression for each 11 (Q^) is finite and it is only the simple expansion technique which 
creates apparent divergences. Our results are 

^n-W) = 4S(-Ti + ...)-^(TiQ^s(Q^) + f + ...U0f7^'l (21) 




""''fl^'iQ') = ^(Ti + ...) + ^(Tig2S(Q') + -^ + ...)+0( — ) (22) 

^•°n'^W) = ^(-^i + ---) + ^(g^j (23) 

^n-(g^) - ^(-3T. + ...)-^(^gW)-T. + ...) + o(i,)(24) 

^U'^\Q') = Oj^i^^ (25) 

We have defined 

M""^ = Ti [m {X" , X^}) (26) 

where |A'^,A''> is an anticommutator. The ellipses at a given order in l/Q"^ denote the 
undetermined contributions from divergent integrals which exist at higher order in the 
naive integrand expansion. (We will discuss the calculation of these terms below.) We 
have not shown the current mass matrix Ai except for the leading logarithmic terms 
at order Q~'^{lnQ^Y, where we have retained the terms linear in A4. The T„, are one- 
dimensional integrals which diverge logarithmically and require renormalization. When 
(|^ is used to evaluate these integrals from the QCD scale A up to a cutoff Q^, we obtain 



47r2yA2 x + S2(x) 3(1 -2rf) ^^^' \ A^ ^ ' 

To leading order in ag, the OPE in QCD gives|I|] 

-n-w^) = ^(5,) + o(ij) (31) 

Sfrabfri'^^ _ 3A^" /7,^\ , _^! /^, r-^ rrr-Mi'V I V4™I . -aft /- ^r^ri-r^i^ 
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We have defined 

^'^''(r^.J = vra, (v^r^...r^A>V^r^-r'-AV) (34) 

B'"' = vra7^7M^^{A^A^}^ J] g7'^A^g\ (35) 

\ q=u,d,s,... I 

Recall that q denotes a single quark fiavor whereas tp is the quark A^j-plet. The r'' are 
SU(3) color generators normalized by Tr (r^'r*) = 25'^'*. 

A comparison of (|21| - p5D and (|29| - p3D gives our expressions for the leading vacuum 



condensates. At the outset of this section we neglected (agG'^^^G^'^'^) contributions to the 
quark propagator, and we now see explicitly that (ugG'^^^G'^^^j is absent from our results 
for the correlators. The results verify, for all correlators except possibly the pseudoscalar, 
that the model correctly generates no dimension two condensates. The integrand of the 
pseudoscalar expression could not be expanded in our simple manner beyond dimension 
two without introducing quadratically divergent integrals. Whether or not the offending 
terms actually sum to zero cannot be determined with this technique. At order 1/Q^, the 
terms in the vector, axialvector and scalar correlators which involve the current quark 



masses also agree with the OPE if we can neglect the ellipses in ( pT| - p^ ). (We will jus- 



tify the neglect of these terms below.) The successful result for the scalar correlator in 
particular is a reflection of the scalar vertex present in the GNC Lagrangian. 

At dimension six, we neglect Ai corrections and find expressions for the 4-quark 
condensates that appear in each of three correlators. 

9^"*(7^) + i3"* = -18Tig2s(Q2)5afe _ 3X2^"^ + . . . 
9^"'(7m75) + S"' = 18T,Q^J:{Q^)5'''' + QT26''' + ... 

Recollection of the Q^ dependence of (qq) given in (||) reveals that these functions ex- 
hibit the known (InQ^)^'^"^ dependence on (5^.[0] However, the presence of undetermined 
contributions at the same order in l/Q"^ prevents us from using this approach to obtain 
complete expressions for these condensates. To proceed, we examine a specific example 
in detail. 



Rigorous 1/Q^ expansion for a specific mass function 

The loop integrations from Fig. |I| can be performed rigorously for the simple mass 
function which has typically been used in low energy GNC calculations. 

mo ~ 300MeV represents the scale of the constituent quark mass and the data require 
2 < y4 < 3.0 Our goal is to build a mass function which becomes (|39|)[([^)] in the limit of 
small[large] Q^, but we begin with a discussion of (B^) alone. 



For any A > 1.44, the Euclidean propagator can be expanded in a convergent geometric 
series. 

1 1 ^ Mm2-S2(g2)\*^ 



g2 + S2(g2 
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(40) 



(41) 



.-o^t^^K^-O! (g^ + Am2)'^+2m 
When this relation is used for each propagator in the expression for a two-point correlator, 
the result contains four summations, but each term in the infinite series is an integral that 
can be evaluated by introducing a Feynman parameter in the standard fashion. The final 
expression does not, in general, resum to a simple functional form, but can be reduced to a 
single summation over hypergeometric functions which can then be summed numerically 
to any desired accuracy. 

For example, the leading nonperturbative behavior of the vector correlator is 



V-Frab 



W\Q' 
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(47r): 



M'Hu 
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-Yv{v 



+ 0{v?) 



where IVff ) denotes the infinite series and 



u 



V = 



A+1 
A^rriQ 
Ami 



g2E(g2 
E(g2^ 



as g^ 



oo 



g2 s(o) 

(A +1)2 S3(0) 



as g^ 



oo 



as g^ 



CXO 



A3 Q^^iQ^) 

For any A > 1.44, Yv{v) is well-approximated by the following simple form. 



1 
2 



(42) 

(43) 
(44) 
(45) 

(46) 



The value of Yv{0) is exact, and all derivatives at v=0 are infinite in the true (series) 
expression as well as in this approximation. 

If we assume (incorrectly) that this result remains valid when (|39D is modified to be 
consistent with the correct asymptotic form of (|^), then we get 



V-Frab 



W{Q' 



4^^°'<'44! 



-1 r 



In 



.^2(0) 



- Yvd) 



+ Oiu' 



V = 



s3(o) ^ -3S3(o) / g2 
g2E(g2) Ad'K^qq) y" K^ 



(47) 
(48) 



A comparison with (^) shows that we have obtained the correct leading logarithm except 
for a missing factor of 1/d. This is to be expected. If we had put the true form of the 
mass function (0) into the integral of Fig. |I|, there would have been an extra factor of 1/d 
which arises from the integration. This is easily demonstrated by computing the large Q^ 
behavior of Ti, as defined in (|27|) , with a cutoff. 



Q , xS(x) 
ax 

A2 X 



dx— 

A2 X 



•4) 



d-l 




(49) 



It is now clear that we were justified in neglecting the uncomputed terms containing Ai 
in ( pll - piD as was postulated in the previous section — the uncomputed terms produce 
functions Yv{v), Ya{v) and Ys{v) which are logarithmically suppressed relative to the 
computed term. 

We now set A^=0 and consider the leading nonperturbative contributions to the vector, 
axialvector and scalar correlators. The pseudoscalar correlator could also be obtained, 
but the presence of the pion propagator (Figs. and ^ makes the calculation rather 
tedious. The pion contribution to the axialvector correlator is slightly more pleasant, and 
the transverse nature of this correlator provides a valuable check on the calculation. The 
relation between the vector minus axialvector correlator and the pion mass difference, 
which we discuss below, is another incentive for pursuing the axialvector calculation. Our 
results are 






346' 



ab 



(47r)2 
~{4^ 
(47r)2 



M^f 
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- Za{v 

- Zs{v\ 



+ 0{u^) 
+ O(m^) 
+ O(n^) 



(50) 
(51) 
(52) 



where Zx{v) denote the infinite series and u, v are defined in ( ^4| -| 
the Zx{v) are well-approximated by the following simple forms. 

67 
Zv{v) ^ — + 0.24v/^ 
51 

Za{v) ^ ^ + 0.40v^ 
5 Jv 



). For any A > 1.44, 

(53) 
(54) 
(55) 
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The values of Zx{0) are exact, and all derivatives at v=0 are infinite in the true (series) 
expressions as well as in these approximations. 

If we consider a modification of (^) to make it consistent with the correct asymp- 
totic form of (^, then we see that each correlator regains the known leading logarithmic 



behavior of (|36| - pq ). Moreover, we can now estimate the size of these condensates. 



9^"'(7^) + S'^' 



28d^7r^S''' {qqf (\n^ 



r,/n2 
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68^2 vr^^"' {qqf I In 
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- Zv{v 

- Zs(v] 



(56) 
(57) 
(58) 



We may compare these to those in the vacuum saturation approximation: 
9A'^\^,) + B^' = ^-^drr'S^' {qqf L^) ' 



9^"''(7m75) + B' 



ab 



^<i.=«-(5,)'(ln| 



2\ -1 



iab 



SA'^'ia,,) + B' 



ab 



fiu J 
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dn'6'^'' {qqf ( In^^^ 



-352 , orab /^„\2 / ,„<5 



A2 



(59) 
(60) 
(61) 



The two expressions have the same leading logarithmic dependence on Q^, but the numer- 
ical factors in front differ. We could make similar remarks here as made in the discussion 
surrounding (^7|-^). 

Of greater interest are the terms which are subleading to the leading logarithm, and 
which correspond to terms which are often neglected in approximations to QCD such as 
vacuum saturation. We may compare the terms in brackets in (^) and ( p^ - |58| ) to the 
leading logarithm, \n{Q'^/A'^). By making reasonable estimates of the various quantities 
involved, 



lOOMeV < A < SOOMeV 

400MeV < S(0) < 500MeV 

-(230MeV)3 < {qq)^(\nl^y'^ < - (180MeV)3 



(62) 
(63) 
(64) 
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and by setting (i=4/9 and Q=lGeV we obtain 



2.4 < In (f^) < 4.6 (65) 



2.2 < ln(^)-lV(t)) <4.2 (66) 

1.7 < \n{^)-ZAiv) <2.3 (67) 

0.6 < \n[^)-Zviv) <1.5 (68) 

-12.1 < In (^) - Z5(t}) < - 10.8 (69) 

The difference between the leading logarithm, ln((5^/A^), and our actual result is small for 
the 1/Q^ term (vector correlator), but larger for the 1/Q^ terms in the various correlators. 
In the case of the scalar correlator the difference is remarkably enormous. The leading 
logarithm still controls the Q^ dependence, but the magnitude of the implied condensates 
at a given value of Q^ has been shifted substantially. We discuss the implications in the 
conclusion. 

The 77'^ — 77^ electromagnetic mass difference 

Of most physical interest is the vector minus axialvector, or left-right (LR) correlator, 
which is a direct measure of chiral symmetry breaking. For 7W=0, 

^^I[{Q^)6''^ = ^n"^(Q2) _ ^fl"''(Q2). (70) 

We now return to our full expression for the LR correlator, without a l/Q"^ expansion, 
and use it to calculate the vr"*" — 7r° electromagnetic mass difference. This mass difference 
is given by@ 

Ami = Pr.r dQ' \Q' ^^n(g')l (71) 



4vr/2 Jo 

where a is the electromagnetic coupling. We will neglect the tiny effects of nonzero current 
quark masses. As discussed earlier, the large(small) Q^ portion of the correlator will be 
obtained from Fig. |l|(Fig. ^ plus meson loops). In fact, it will become evident as we 
proceed that we can get an upper bound on the mass difference by simply using Fig. |I] 
for all Q\ 

We must also decide how to extend the asymptotic form of S((5^) given in (0) to 
smaller Q^. The following simple ansatz contains four parameters: Atuq, Bitiq, C and 

12 



M^. (As in 01, mo represents the scale of the constituent quark mass. Since it is not an 
independent parameter here, we are free to choose it to be numerically identical to its 
value in 1^.) 



Bml 



[Ami + Q' 



l + C{ln{l + ^)y~ 



(72) 



The motivation for this functional form comes from the M —>■ oo limit, where ^(Q^) 
becomes the mass function used originally in the low energy GNC model. We wish to 
include the correct logarithmic behavior of the mass function at high energies without 
making significant changes to the original low energy form. The ten dimensionless quan- 
tities Li{fi) which appear in the standard chiral Lagrangian]^ have been obtained from 
the GNC model (i.e. Fig. ^ in the M — > cxd limit. 01^ They are sensitive to the shape of 
the mass function, but not to the overall scale, so from this we can determine the value 



of Ami- 



(470 MeVf < Ami ~ (550 MeV)^ 



(73) 



One constraint on the three remaining parameters in ( [7^ ) comes from demanding that 
S((5^) satisfies the high energy behavior of (|^), where the numerical value of (qq) is 
known phenomenologically. This gives 



Bml 



-Mn'^ 



{(l<l)f. 



In 



'A2 



(74) 



Two more constraints arise from the low energy behavior of the LR correlator, which 
must reproduce the experimentally determined coefficients of the chiral Lagrangian. In 
the notation of [0, the relevant coefficients are Fq and Liq. 



Q^^^n(Q^)=Fo^ + 4Q%o(/i) 



6(47r): 



In 






+ 0[Q') 



(75) 



Nf is the number of quark flavors. The dependence of Liq on the renormalization scale /i 
is canceled by the logarithm (which comes from internal meson loops), so that ^^Tl{Q^) is 
independent of yU. We will eliminate the parameters Bm^ and C by using the phenomeno- 
logical values for (gg) and Fq. The final parameter M will be determined by requiring 
the mass of ([72|) in the calculation of Fig. |^ to produce a value for Liq which is within, 
say, 25% of the original GNC result. 
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Before proceeding, we point out an interesting relation between the high and low 
energy calculations. From the evaluation of Fig. |^, we obtain 

and a more lengthy expression for Liq. It turns out that exactly the same expression for 
Fq happens to come from Fig. |l]. The two expressions for Lio are not the same; Fig. |l] 
correctly predicts that Liq is negative, but the magnitude is only about half of the correct 
GNC result. This means that if Fig. |I| is used for all momenta, the slope of the LR 
correlator at Q^=0 is too shallow (see Fig. |^), and an upper bound is obtained for Am^. 
We choose the strongest upper bound by using the smallest value of M which keeps Liq 
(of Fig. D within 25% of the original (M -^ oo) GNC result. 
For definiteness, we use (i=4/9 with the following inputs, 

Fq ^ 88MeV (77) 

(gg)^(^ln^^ ^ -(220MeV)=^ (78) 

Note that this last expression is independent of /x. With these values, the parameters B 
and C of ( [7^ ) are of order unity, and the scale M is about 3 GeV. Fig. |T] then gives an 



upper bound on the electromagnetic pion mass difference. 

m^+ -m^o, < 5.1 MeV (79) 

To obtain a direct estimate we will calculate the high energy expression for the LR 
correlator from Fig. |l] down to some intermediate scale Qligh and use Fig. |] plus meson 
loops below the scale Qfow such that 

Q' "^n(gL) = Q' ^^n(QL,.) = R' (so) 

Between Qf^w ^^^ Q\ighi ^^is function will be approximated by the constant R^. The LR 
correlator is plotted versus Q^ in Fig. ^, and our result for the pion mass difference is 
shown in Fig. ^ as a function of R. The experimental value (after subtracting md — rriu 
effects) isg 

[m,+ - m^ol^p^ = 4.43 ± 0.03 MeV (81) 
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and corresponds to 50MeV < R < 55MeV. This implies that for this calculation, our low 
energy model is good up to a scale of order Qf^^ ~ 400 — 450MeV and our high energy 
model is good down to a scale of order Qlig^ ~ 750 — 850MeV. Both of these scales are 
very reasonable. 

Conclusions 

One consequence of the dynamical breakdown of chiral symmetry in QCD is the gener- 
ation of a momentum dependent light quark mass. In the context of the vector, axialvector 
and scalar two-point correlators we have shown how this effective mass can be included 
systematically in calculations. The Ward-Takahashi identities and chiral symmetry are 
respected, as are the first few terms (at least) of the operator product expansion (OPE). 
An interesting result is the existence of terms at order 1/Q^ which do not contain the 
leading logarithm but which are not insignificant, especially in the case of the scalar 
correlator. 

It appears that for the scalar two-point function a naive application of the OPE in con- 
junction with the vacuum saturation approximation does not adequately describe some 
expected physics of QCD, namely the physics associated with the dynamical quark mass. 
This is perhaps not surprising. Practitioners of sum rules have long claimed that there 
is something deficient in the usual application of the OPE in the case of the scalar and 
pseudoscalar two-point functions. This is particularly clear in the pseudoscalar case, for 
which the conventional OPE does not adequately account for the pion appearing in the 
sum rule. This has led to speculations of additional contributions to these OPEs, includ- 
ing instantons0, renormalons ||1 0|| , and effective four-fermion interactions 0]. Our work 



indicates that the additional contributions will also have to reflect effects associated with 
chiral symmetry breaking, and in particular the momentum dependence of the dynamical 
quark mass. On the other hand in the case of the vector and axialvector two-point func- 
tions, the subleading terms in our analysis are not substantial enough to say that there 
is a serious conflict with standard treatments. 

We have described in this work a minimal model. It could be extended for example by 
including an effective wavefunction renormalization parameter in the quark propagator, or 
by adding extra terms to the vertices which maintain the Ward-Takahashi identities and 
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the OPE. Interestingly enough, we find that the minimal model is sufficient to account 
for the vr"*" — 7r° electromagnetic mass difference. 
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Figure 1: The complete QCD contribution to a two-point correlator. Both propagators 
are full propagators; one vertex is full and the other vertex is a bare ^[75] or ^7^ [75]. 



Figure 2: The quark contribution to a two-point correlator in the GNC model. Meson 
contributions are not shown. 



Figure 3: The two distinct components of the full axialvector and pseudoscalar vertices 
as derived from the GNC model. The dashed line represents a pseudoscalar meson prop- 
agator which is generated from the diagrams of Fig. 11. 



Figure 4: The GNC Lagrangian does not contain an explicit meson propagator, but a 
propagator is generated by the quark loop diagrams shown here. 



Figure 5: The vector minus axialvector two-point correlator obtained from the GNC 
model at low energies and from Fig. |^ at high energies for a typical choice of parameters. 
The two pieces are matched to a constant, i?^, in the intermediate region. 



Figure 6: The electromagnetic mass difference of the pion as a function of R, defined by 
(pop, for a typical choice of parameters. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411412v2 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411412v2 
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